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Abstract 

A new method to find the propagation equation system governing the scattering 
of an electromagnetic wave by a nonlinear medium is proposed. The aim is to let 
the effects appear spontaneously, deleting as far as possible the phenomenological 
ideas. In this way, we obtain propagation equation systems that encodes several 
nonlinear effects. Once these systems obtained, the numerical values of the tensors 
that characterize the answer of the medium to an electromagnetic perturbation give 
weights to the different effects. 

This aim is partly reached in this study, especially when treating harmonic gen- 
erations. For this, we start from the Maxwell's equation and give rigorously all 
the hypothesis needed to attain equation systems that can be solved, at least from 
a numerical point of view. Finally, a symmetry on the susceptibility tensors that 
ensures that a medium is lossless from the electromagnetic point of view is worked 
out. 

1 Introduction 

The usual route to treat nonlinear media from the electromagnetic point of view is to 
expand the electric permittivity as a function of the electric field. These expansions 
quickly get huge and so one has first to select the effect we want to consider, and then 
systematically simplify the equations, keeping only the terms that contribute to this 
effect ((H |2])- A look at the literature shows that the precision that can be obtained 
with this method is out of question. 

Nevertheless, it would appear more satisfactory to start only with Maxwell's equa- 
tions and to let the effects appear by themselves, introducing a systematical way to 
simplify the equations. The simplification is thus a priori, generic ; the relative im- 
portance of the observed effects being only a posteriori determined, by the numerical 
values of the susceptibility tensors. This paper aims at drawing this new road. 
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For this, some general work on the constitutive relations had to be done. We ex- 
pose it in the subsection 12.11 Then, with some assumptions, we recover in the next 
subsection the standard expression of the electric susceptibility tensors. Most of these 
assumptions are justified for the presentation of this paper, and avoiding them would 
not have lead to insurmountable difficulties. 

Having expressed the answer of a medium, we obtain the propagation equation 
system satisfied by the electric field in l3.ll This system, parameterized by a continuous 
parameter, is far too complicated to be solved, even with numerical methods. This 
leads us to partly leave our general aim, and concentrate on harmonic generation. In 
this context, we introduce a new notion, called the degree: it determines, in generic 
terms, which interaction between the Fourier components of the electric field have to be 
taken into account. The subsections l3.3l and l3.4l are devoted to the propagation equation 
systems in the lowest order of nonlinearity and in the lowest degree; the physical effects 
are then clearly identified. We hope the obtention of these equation systems from ab 
initio principles will convince the reader about the interest of the method. 

In a second part (section 0J, we try to answer, we think more rigorously that what 
can be found in the literature, to the question: " how the electric energy variation is 
determined by the susceptibility tensor fields of a medium?' If the general answer is 
still unknown, we give a sufficient condition for a medium to be lossless. 

This paper is intended to be appreciated by the theoretical physicists, either work- 
ing in the nonlinear optics field or not, as well as the applied physicists. To this aim, 
we start from general concepts and end with equation systems that can be numeri- 
cally solved. Having clarified the way these systems appear, we have realized some 
simulations, one of them appearing in the companion paper Q- In particular, these 
simulations confirmed our prediction on the energy criterion. 

2 The electromagnetic constitutive relations 

2.1 A tractable expression for the inductive fields 

As said in the introduction, the aim of this article is to expose a new route to the equa- 
tions of nonlinear optics, starting only with the Maxwell's equations. We will work at 
a completely classical level, even if, of course, it is hoped that the effective characteris- 
tics of the medium will in future works clearly come from a microscopic point of view. 
In other words, in this study, the characterization of a medium is considered to be fully 
determined as soon as a set of (susceptibility) tensors are given. 

Let us start with the Maxwell's equations. We consider here that the primitive 
fields are the electric field E and the magnetic field B. The first set of equations gives 
a 'coherence' condition between E and B: 



The second set of equations links the induction fields with the sources of the electro- 
magnetic fields. We thus denote by D the electric induction, B the magnetic induction, 




(1) 
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p the charge density and J the current density. Then, the second set of Maxwell's 
equations is 



( 



V x H <9 t D 
V- D 
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Supposing that the sources are known, we are looking for the electromagnetic 
fields. It appears that we have to add some hypothesis to answer this problem, for, 
in each point s of K 3 and for any time t, we have twelve unknowns (the three com- 
ponents of the four vectors E(s, t), B(s,<), D(s, £) and H(s,i)), but the Maxwell's 
equations give only two vectorial and two scalar equalities. The set of Maxwell's equa- 
tions being independent of the medium in which the fields oscillate, we have to encode 
the (electromagnetic) characteristics of the medium, with two (vectorial) relations that 
give the inductive fields in function of the primitive fields. We pose the following set 
of equations: 



R respectively. As said above, the 'fabrication' of the electromagnetic induction is done 
at another scale than the one considered here, and we are reduced to see it like a black 
box: at the input are, for each place s and each time t, the primitive fields, and s and 
t; at the output are the inductive vectors D(s, t) and H(s, t). Several authors (0[5), 
see also [6|) give expressions similar to (fJHUl, but we think that the one we present is 
more rigorous and allows a clear transcription of the properties of the medium. The 
article [7|, and some references found here, is not restricted to pure electromagnetic 
phenomena. 

The transcriptions, between the properties of the medium and the restriction it in- 
duces on the functionals 35 and Sj are the following one: 

• The causality of the medium implies that 



which means that the times r at which are evaluated the electric and magnetic 
fields that contribute to D(s. t) and H(s, t) are prior to t. 

• The locality (in space) of the medium implies that 



Allowing the notations' abuse that consists in keeping the same names for the 
inductive functionals, the new expression are 




a 



(i) C (-00, t], (a, b) G {e, b} x {d, h} 



S b a (s) = {s}, (a,b)e{e,b}x{d,h}. 
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f U(s,t) = D({E( S ,T)} TeT 4 {t) ,{B( S ,T)} TeT , {t) ,s,t) 
\ H(s, t) = S)({E(s, T)} reT H (t) , {B(s, r)} reT h (t) , s, t) . 



• A medium is nonbianisotropic if 



<9 B S = 
d E Sj = 0. 



In this case, we have 



D(s,i) = 3)({E(o-,r)} (0 . )T)eS |(s)xTd(i),s,t) 
H(s,i) =i5({B(<T,T)} (<y . iT)eS h (s)xT h (t) ,s,i) 

• A medium is homogeneous in time, or, as is commonly said, stationary, if the 
functionals do not explicitly depend on the time t: 



f dt® = 
I d t Sj = 0. 

The constitutive relations then get: 

D(s, t)=D ({E(cr , r)}( £r)T ) e 5d( s ) xT d( t ) , {B(cr, t)} ((T , r ) e s^(s) xT 6 d (t) > s ) 

H(s, t) = #({E(<T , r)}( - )T ) e 5h( s ) xT h( t ), {B(<T , T )}( C r,r)eS 6 ft ( S )xT b ''(t)> s ) ■ 

We would like to stress that this hypothesis is, in a way, delicate, when treating 
nonlinear optics. Indeed, for the nonlinear effects to be important, the intensity of 
light has to be large. But if this intensity is too large, then we cannot neglect the 
effects of the electromagnetic field on the medium. Some effects like saturation, 
damage, etc. do appear, and we leave the realm of stationary media. 

The transcriptions of locality in time (commonly called nondispersive) or homo- 
geneity in space in the inductive functions should be obvious. We can then define 
the (local) electromagnetic vacuum: a vacuum holds at a point s and a time t if the 
following equalities of vectors are satisfied: 



( D(s,t) = e E(s,t) 
I H(s,t) = l/MoB(s,i). 

This leads to define the 'answer' of a medium to an electromagnetic perturbation as the 
difference between the induction in that medium and the one in a vacuum: 
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f P e :=D-e E 
\ P m :=H-l//ioB. 

These two vector fields are called the electric polarization vector field and the magnetic 
polarization vector field - we prefer to write P m for what is usually denoted by M 
(termed the magnetization vector field) for the coherence of this article. We now arrive 
at a key step. In the following definition, we denote by V the dual of the vector space 
V. 

Definition 1 Smooth Medium 

A nonbianisotropic medium is called smooth in a neighborhood of a point s and a 
time t if the electric and magnetic polarization vector fields admit Taylor expansions 
of the following kind: 

f P e = P<°> + ]T Pi") 

nGN 

p _ p(0) , V p(n) 

with 

s p OO />OG 

P e n) M) =£0 / dti--- dt n Q< n >(s,t;si,ti,--- ,s n ,t n ) 

J— oo J —oo 

E(si,ii)---E(s„,i„) 

S /" OC /-OO 

P^ ) (s,t)=£o/ (fti---/ tft n Q^ ) (s,*;si,ti,--- ,s n ,i n ) 

J —oo J — oo 

H(si,ii)---H(s„,f„) 

Q(") : M 4 x R 4 ™ M 3 ® (R 3 )'®" 
: M 4 x R 4 ™ -> M 3 (R 3 )'®". 

A similar definition can be applied to bianisotropic media as well: the tensors Qe 
and Qm' have then to be contracted with the 2™ combinations of E and B. Though 
treating these media, among which we find the chiral ones, do not bring important new 
difficulties, we won't anymore be concerned with them, to keep formulae of 'reason- 
able' size. 

Due to the nonbianisotropy hypothesis, the electric and magnetic inductions can be 
treated separately. We will from now on concentrate on the electric one. Hence, when 
no confusion is possible, we will write P in place of P e , is place of Qi n \ say 
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the polarization vector instead of the electric polarization vector, etc. The treatment 
required for the magnetic part is completely similar. 

The set of operators Q^ n \ when applied to (s, t; s\, t%, • • • , s n , t n ), give the ef- 
fect, at the point s and time t, that the electric field, evaluated at the points s% , ■ ■ ■ , s„ 
and times t\, ■ ■ ■ ,t n have on the electric polarization vector at the point s and time t. 
The zero-th order of the electric polarization vector PW(s, t) corresponds to a spon- 
taneous nonzero electric moment. It does not vanish for ferroelectric materials. For 
convenience, we write p(°)(s,t) = eoQ^°\s,t). The term P'"', for n S N is called 
the n-th order of the polarization vector. 



2.2 A tractable expression for the polarization vector 
2.2.1 The definition of the electric susceptibility tensors 

The expansion of the (electric) polarization vector given in the definition (Q]) is still too 
complicated, and needs further simplification to be studied. The following hypothesis, 
still motivated for the presentation, is that the medium is spatially local. The electric 
inductive functional has the form 

D(s,t) = 2)({E(8,r)} T63?(t)> 8 I t) 

so that 

dh--J dt n Q^ n \s,t;t u --- ,t n ) 

E(s,ii)---E(s,t„). (5) 

Once again, we have adapted the functionals without changing the name (now, we have 
Q( n ) : R 4 x K" M 3 ® (K 3 )'®"). 

Of more fundamental nature, let us suppose that the medium is stationary; this 
means that, first, the spontaneous polarization does not depend on time, and secondly, 
that for any point s and for any duration time T, the effect of the electric vector field 
evaluated at the point s and at the times t\, ■ ■ ■ ,t n on the polarization vector field 
evaluated at the point s and at the time t is the same than the effect of the electric vector 
field evaluated at the point s and at the times t\ — T, • • • , t n — T on the polarization 
vector field evaluated at the point s and at the time t — T, i.e. 

QW( s ,t)=Q(°\s,t-T) 

and 

QW(B,i;ti,..- ,t n ) = Q( n \s,t-T;h-T,--- ,t n -T), Vn G N, VT e R. 
When considering stationary media, a usual trick (|2l) is to define 

: M 3 -> M 3 , 

s h+ Q(°)(s,0), 
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i?W : l 3 xM" -> R 3 ® (R 3 )'®", 

(s, ,i n ) h+ Q( n )(s,0;-ti,--- ,-i„), neN. 

This function i?( n ' is called the response function of the n-th order, because it allows 
to write the relation (0 as a convolution on R™ between R^ and 

E (,l) : R 3 x W -> (R 3 )®", (s,ti, • ■ • , t„) i-4 E(s, h) ■ ■ ■ E(s, t n ) : 
pW(s,t)=e o (i2 (n) (8,-)*„EW(8,0)(*) 

/oo /*oo 
tfti-- - / dt„flW(s, ,t-i n )E(s,*i)---E(s,t n ) 
-co •/ —CO 

In the last equation, we make the transformation t j \-± t — tj to obtain 

/OO /» oo 

dti--- / dtnR^is,^,--- ,t n )E(s,t-ii)-.-E(s,t-t n ) 
-oo J — oo 

and then we express the electric field in its Fourier basis, 

/CO 
dw i e-*-^ t - t ^E(8 J w i ), J e {!,-■■ ,n}. 
-co 



This leads to 



/CO />C 
dwi e-' 1 " 1 * • • • / 
-OO J — ! 



X (n) (s,wi, • • • ,cj„)E(s,wi) • ••E(s,o; n ) 

/OO p OO 

-OO — OO 

E(s, Wl )---E(s, W „)e-^ 1 +-+"") t (6) 

where 

X (n) :R 3 xR" C 3 ®(C 3 )'®", (7) 

/OO /*CO 
dtf- / dt n R( n \s,h,--- ,t„)e iUM 
-oo J —oo 

with a; = (wi , ■ ■ • , w„), t = (ti, • • • , t n ) and the usual duality product on R™. We see 
that x^> which is a tensor field of rank n + 1, called the (electric) susceptibility tensor 
of the n-th order, is the Fourier transform (up to a scalar multiplication by (27r) Tl ) of 
the response function Ry 1 ' : 



X< n >(s, Wl) -.- ,u; n ) = (2n) n R^(s,uj 1 ,--- ,u n ). 

X^(s,ui, ■ ■ • ,uj n ) gives the way the components, evaluated at the point s, of 
the electric field oscillating with the angular frequencies u)\, . .., u) n contribute to the 

polarization vector oscillating at the angular frequency cji H \- CJ n at the same place 

s. The unit of x {n) is V 1 ' 71 . 

By symmetry, we define x^ by 

X (0) : R 3 -> M 3 , 

s i-> i? (0) (s). 



2.2.2 The symmetries of the susceptibihty tensors 

The susceptibility tensor fields exhibit several symmetries. These ones are important 
to better understand these tensors, and will allow to greatly simplify the results of the 
following sections, about propagation equation systems and energy criteria. 



The intrinsic permutation symmetry We will now present the intrinsic permutation 
symmetry of the susceptibility tensors. Although the result of this paragraph is well- 
known (|[I||2][8)), we demonstrate it because we think that the way we derive it is more 
systematic. For this, let us develop the equation (|6]l in components: 



/oo />oo 
-oo J — oo n 



—i{ui\-\ \-uJ n )t 



where the summation convention over repeated indices is used. The intrinsic permuta- 
tion symmetry is based on the fact that the x 1 ■ are not uniquely defined: for 
j in {1, • • • , n}, the indexes ij and the variables Wj are dummy, so for any bijection a 
from {1, • • • , n) to itself, we have: 

/oo />oo 
diJal--- / du an X {n)l . . (B,W d ,.-. ,U m ) 
I *crl ' ' ' Ian 

-oo J —oo 



E i " 1 (s,u <yl )---E l " n {s : uj an )e 



hi*V„)t 



Thus, denoting by S n the symmetric group on a set of cardinality n (we recall that this 
group is of order n!), we have 

poo poo 

E 1 ' 1 (s, u al ) ■ ■ ■ £>™ (s, w CT „)e- l( "" 1+ - + ^" )t . 
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But, loosely speaking, duj al ■ ■ • du an = dux • • • duj n and 



-iujjt 



for any a in S n . We thus have 



/OO />00 i 

-°° J -°° <xeS„ 

£ 41 (s, wi) • • • E l " (s, w„)e- 4( " 1+ -+"" ) *. 

This leads to define the (not underlined) tensor \^ — e i X^ % i ± ... i n ® en ® 1 " 1 ® 
e 1 " , once again called susceptibility tensor of order n, by 

x (n) : I 3 x M™ -> C 3 ® (C 3 )'®", (8) 
(s,wi,--- ,w n ) i ^ V e;x ( " )4 , , (s,cj ct i,--- ,w CT „) 

H< ' * <-crl "' tan 



and to deduce the following expression of the polarization vector of the n-th order: 

/oo />oo 
dux-- - / dw„x (n) (s,wi,-" ,w n ) (9) 

E(s, Wl )...E(s, W „)e-^ 1 +- + "") t , 
with a symmetric in the sense that 

i Tl ...» Tn (s,w r i, • • • ,w Tn ) = x ( " )j ^...^(s.wi, • • • ,a; n ), Vt e «S n . 

This symmetry of the (not underlined) susceptibility tensor is called the intrinsic 
permutation symmetry. It allows for example to consider only sets of angular frequen- 
cies (u>i, • • • ,L0 n ) ordered such that uj\ ^ u; 2 ^ • • • ^ oj n . 

Finally, since the zero-th order susceptibility tensor does not present degeneracy, 
we define '■— X^ ■ 

The Hermitian symmetry The electric field (and the polarization vector field) E 

being real, the harmonic components satisfy E(s, to) = E(s, -co) for all s in R 3 . This 
leads to the Hermitian symmetry of the susceptibility tensor fields: 



X^(s,cj!,--- ,w„) = x (n) (s,-u} 1 ,--- ,-uj n ). 
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3 The propagation equation systems 



3.1 The general propagation equation system 

The magnetic response The preceding section shows how the electric answer of a 
medium is encoded in the (electric) susceptibility tensors. A similar treatment shows 
how the magnetic answer of a medium is encoded in the magnetic susceptibility ten- 
sors. Nevertheless, for the sake of simplicity of this paper, and also because at the 
wavelength we are looking at, nonlinear magnetic effects are usually negligible, we 
suppose from now on that the magnetic characteristics of the media are smooth, linear, 
local in space and stationary. This means that the magnetic inductive functional has the 
form 

H(s,f) = #({B(s,T)} T6T h (t) ,s) 

so that 

/oo 
cL;^ ) (s,w)B(s,a;)e-- i . 
-oo 

To recover the usual conventions, we define the relative permittivity by: 

£r(s,w) = (l + X^M)" 1 , 

so that 

B = /xH. (10) 

This last equation closes the sets of Maxwell's equations. Reporting the expres- 
sions of the inductive vector fields in terms of the primitive fields ((0 as well as the 
development of the polarization vector field in the definition[TJfor the electric response, 
( [Tol l for the magnetic response) in the first set of Maxwell's equations leads to the 
propagation equations. This section is devoted to present them, and to suggest some 
approximations in order to be able to implement them in a computer program. 

The general propagation equation system On physical grounds, we suppose that 
the permeability tensor never vanishes so that the two vectorial Maxwell's equations 
lead to 

V x (/2 _1 (s,w)V x E(s,w)) - oj 2 I)(s,uj) = iwJ(s,u) (11) 

for all lo in R and s in R 3 . 

The dependence of D in E and in the (electric) characteristic of the medium were 
obtained in the first section, so that we just have to insert them in this expression. For 
this we have to write the polarization vector in the Fourier basis: 
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1 f°° 

pW(s,w) = — / dtPW(s,t)e 

27r ./-oo 

/•oo 



,■00 



2tt 

E(s,cji) • ••E(s,a; n )e" 

/■ OO 

du-L--- I dw„x (n) (s,wi, • • • ,u n ) 



so 



— OO J —oo 

(loi H hw„ — u;)t 

oo /*oo 



— oo 



E(s, wi) • • • E(s, w„)<5(wi H h w n - w) 

where <5 is the Dirac distribution; the last equation follows upon integration upon the 
time variable. For the spontaneous polarization vector, we have 

p(°)(s, W ) = e X (0) (s)^H- 
The propagation equation system is therefore 

V x (/i- 1 (s,w)V x E(s,w)) - w 2 £ (e(s,w) + x (0) (s)S(lj) (12) 

/oo /*oo 
du>i ■■■ / dw n x (n) ( s > w i) ' • ' ,w„)E(s,wi) • • -E(s,w„) 

(5(^i + ■ • • + ui n — ujfj = zu;J(s, (J). 

3.2 Nonlinearity of the second order 

The introduction of some notations Of course, the propagation equation system 
(fT2l is too complicated to be solved directly. We will thus suppose from now on that 
no nonlinearity higher than the quadratic one are present. We will treat this order of 
nonlinearity with some detail; the other ones being direct but more and more sophis- 
ticated generalizations of this case. Also, though its presence does not lead to high 
difficulties, we will neglect the spontaneous polarization. The propagation equation 
system thus becomes 

V x (/i~ 1 (s,w)V x E(s,w)) 
-w 2 £ (e(s,w) + / dw 1 x (1) (s,wi)E(s,u 1 )^(wi-w) 

J — oo 

/OO POO 
dujx / dw 2 X (2) ( s ! w i! w 2)E(s,wi)E(s,w 2 )<S(wi + w 2 - w) 
-oo J —oo 

= zwJ(s, U)) 

for all u in K and s in R a . 
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In order to facilitate the reading of these equations, we introduce some notations. 
First, the usual linear permittivity: ^^(s, u) := 1 + X^H S ; UJ )- Using this and inte- 
grating the Dirac distributions, we have 

V x (/i -1 (s,a;)V x E(s,w)) - w 2 e (e^ (1) (s,w)E(s,w) 
+ / duix x^ 2 (s,wi, oj — wi)E(s,o;i)E(s, oj — oj\) J = iojJ(s,oj) 

J — oo 

for all oj in R and s in R 3 . 

Then, we define the linear Maxwell's operator by 

TW^E^u;)) := -eo X V x (/T^s.uOV x E(s, oj)) + oj 2 t r {1 \s, w)E(s, w); 
this operator is defined such that, in a linear media, we have 

^ ) (E(s jW )) = ^J(s jW ) 
for all w in R and s in R 3 . Finally, we will use the notation 



|E(s,cji), • • • ,E(s,w n )J := \ (n) (s, wi, • • • , w„)E(s, wi) • • • E(s, oj n )\ 

this term is the contribution to the polarization vector P(s, oj\ + ■ ■ ■ + oj n ) of the 
interaction between the vectors E(s, oj\), ■ ■ ■ , E(s,w n ). The propagation equation 
system is now 



^)(E(s,^))H 
for all oj in R and s in 



dm LE(s,wi),E(s,w - oji)\ = J(s,w) (13) 



The harmonic assumption and the definition of the degree We focus on the case 
where the incident vector field is monochromatic: E J (s,i) = E l 1 (s)e _lWjt +E!_ 1 (s)e lWj 
(we note E*(s) for E l (s,pwj)). To simplify the system ( fT3l . which is described by a 
continuous parameter, we suppose that the susceptibility tensor field is nonzero only 
when it is evaluated on the harmonics of the incident angular frequency 

X W(s, Wl ,-, Wn ) = 

if there is a j in {1, • • • , n} such that no integer p satisfies ujj = pojj. This is an impor- 
tant hypothesis, which implies that we will concentrate on the harmonic generations 
and neglect the Raman or Brilloin scatterings, or subharmonic generations. On one 
hand, we have to admit that this is a serious stretch to the general aim stated in the in- 
troduction. On the other hand, the set of harmonic generations is sufficiently important 
and large to pursue this study. 
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The propagation equation of the p-th harmonic is now (we stop to explicitly write 
the point s where the fields are evaluated - note that no confusion can appear, because 
of the locality of the considered media; this is just a return to equalities between vector 
fields on R 3 and not simply vectors in M 3 - we also write A4 l ^ n for M 1 ? 11 A 

M l ™ (E p ) + (pojj) 2 ]T [Eg, E P - q \ = J P) (14) 

for all p in Z. Since the sources oscillates only at the angular frequency uij (and of 
course —lji), we have J p = J p <Ji p m, We note that the harmonic assumption turns 
the propagation equation system from a system described by a parameter ranging on a 
continuous set to a system described by a parameter ranging on a discrete set. 

Nevertheless, (TBI is still too complicated to be solved: it still contains an infinite 
number of equations. To this aim, we introduce the notion of the degree, whose effect 
will be to obtain propagation equation systems described by a parameter ranging on 
a finite set. No absolute definition of the degree does exist; indeed, the way we will 
simplify the system ( TPfl ) will depend on our purpose. To illustrate the choice we have 
to settle on, we give two definitions, and will argue that these two definitions present 
interesting, but different points of view. 

Definition 2 the degree d\ of the approximation at the n-th order of a monochromatic 
field is 

d x :=min{||(pi,--- ,p n )\\uz) > d x (n) 0;Pi w /, ' • ' ,Pn^i) = 0,VsgM 3 } 

Definition 3 the degree d^ of the approximation at the n-th order of a monochromatic 
field is 

doo := min{||(pi, • • ■ ,p n ,pi H hp„)||/ co (z) > d 

=*- x {n \s;PiUi, ■ ■ ■ ,PnUl) = 0,Vs e M 3 } 

With words, this means that, since we can consider only a finite number of compo- 
nents of the electric field, we suppose that the susceptibility tensors vanish if they are 
evaluated at high enough frequencies. In the di-th degree d, all the terms of the kind 
|_E Pl , • • ■ ,E Pi J vanish as soon as |pi| + ••• + \Pn\ is strictly greater than d. In the 
doo-th degree d, all the terms of the kind [E Pl , ■ ■ ■ , E Pn J vanish as soon as \pi |, • ■ • , 
\p n | or \pi + ■ ■ ■ + p n I is strictly greater than d. In both cases, the electric vector can 
be written, for all s in M 3 and all t in M as 

E(s,i)= ]T E p (s) e -^ 4 

pGZ 
\p\<d 

= E + 2Ke{ Y, E p (s)e-^ t }, 
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since E_ p = E p . 

We hope that these definitions of the degree will become intuitive when reading the 
next subsection. 

3.3 The propagation equation systems in the lowest degree 

We add a new hypothesis: the static component of the electric field Eo vanishes. Keep- 
ing it would not have brought huge difficulties, at least in this theoretical work. We 
note that this prevents us from studying the Pockels effect or the optical rectification. 

3.3.1 In the degree 1 (linear case) 

The cases d\ = 1 and doo = 1 are both: 

This means that, whatever we use the d\ or the d^ degree, the coarser approximation 
that we can do in nonlinear optics reduces to the linear case. 

Remark 1 The symmetry of the propagation equation system 

Due to the Hermitian symmetry of the Fourier components of the electric field 
(E_ p = E p J, the one of the susceptibility tensor (see the paragraph The Hermitian 
symmetry, page\9§ and the one of the permeability tensor ( fi(s, —oS) = /i(s, w) ), it can 
be shown that the system \14\ is coherent in the following sense: i/E p is a solution of 
the propagation equation for E p , then E p is a solution of the propagation equation for 
E_ p . This allows to give only the equations for positive frequencies. 

3.3.2 In the degree 2 (second harmonic generation) 
With the d\ definition In the d\ = 2 case, we have 

M m (Ei) = =^Ji, (15a) 

M l 2 m (E 2 ) + (2w/) 2 LEi,EiJ = 0. (15b) 

Indeed, in this degree d\ — 2, can be contracted only with the components E p 
and E 9 such that \p\ + \q\ ^ 2; since Eo = 0, we also have p ^ and q ^ 0. Hence, 
the nonzero susceptibility tensors are x^ 2 )(— w/, — uf), x' 2 ' (w/, — u>i), X (~ u ii 10 1) 
et x^H^i: ^i). But the second and the third ones do not appear in the propagation 
equation system since we suppose that Eo vanishes, and the first one describes the 
generation of the component E_2, that we directly obtain from the equation that satis- 
fies E 2 . 

This propagation equation system has the advantage that it can be solved treating 
only linear equations: one first solves the equation that Ei satisfies and then the one 
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that E2 satisfies, treating Ei as a source. But this does not mean that the system is 
linear: if the incident field is multiplied by a constant factor m, Ei is scaled by m 
whereas E2 is scaled by to 2 . 

Since the equation that satisfies Ei is linear, it is commonly said that we are in the 
framework of the nondepletion of the pump beam. Indeed, Ei can be seen as a tank: 
the fact that it generates E2 does not change its energy. 

With the doo definition In the = 2 case, we have 

Mi <n (E!)+2u;|LE-i,E 2 J = ^^-Ji, (16a) 

A4 in (E 2 ) + (2w/) 2 LEi,E 1 J = 0. (16b) 

Considering this system, the definition of the degree d^ can appear to be nonphys- 
ical: in the generic case, that is without specifying a material, there are no reason for 
keeping the term |_E_i,E2j when neglecting LEi,E 2 J that could generate the third 
harmonic (we recall that ||E_i|| = ||Ei||). This remark argues in the sense that the 
degree d\ is the most natural one from the physical point of view. Nevertheless, in 
practical cases, the degree d^ can be important since it allows to consider a lossless 
medium - a fact, as seen above, impossible with the d\ definition. The details are given 
in the sectionH] 

3.3.3 In the degree 3 (third harmonic generation by a cascade effect) 

A small number of interactions was considered in the degree 2. In the degree 3, we 
consider the interactions that result from a cascade effect. With the d\ -degree, the 
interactions between E_2 and E2 on one side, and the components E_i and Ei on the 
other side are taken into account. With the doo-degree, we consider all the interactions 
where E p occur for p in {—3, —2, —1, 1, 2, 3}. 

With the d\ definition In the di = 3 case, we have 

X'r(E 1 )+ W/ 2 (LE_ 1 ,E 2 J + LE2,E_!j) = ^^Ji, 

M m (E 2 ) + (2^ J ) 2 LE 1 ,E 1 J =0, 

M l t(E 3 ) + (3 W/ ) 2 (le 1 ,E 2 J + LE 2 ,E!j) = 0. 
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Using the intrinsic permutation symmetry, the propagation equation system gets 

M l i n (Ei) + 2uj [E_i,E 2 J = ^-Ji, 

X! ! " l (E 2 ) + (2 W/ ) 2 LE 1 ,E 1 J =0, 
M m (E 3 ) + 2(3w / ) 2 LE l! E 2 J =0. 



With the g?oo definition In the = 3 case, we have 

A^f i (E 1 ) + 2 W |(LE_ 2 ,E 3 J + |E_i,E 2 j) = ^-J 1; 

A^ 2 m (E 2 ) + (2w / ) 2 (2LE_i,E 3 J + lEi.Exj) =0, 

X 3 m (E 3 ) + 2(3^ / ) 2 LE 1 ,E 2 J =0. 

It may seem curious that we consider the third harmonic generation without taking 
into account the third order polarization. In fact, we advise to the reader not to use 
these systems - they were exposed only to give better intuition on the notion of degree. 

In the second order of nonlinearity, the relevant systems are the ones obtain in the 
degree 2. It is a general fact that the n-th harmonic generation is to be studied in the 
n-th order of nonlinearity and in the n-th degree (at least with the two definition of the 
degree presented here). 

3.4 Nonlinearity of the third order 

The remark given in the second order of nonlinearity and in the third degree argues that 
we have to take into account the third order polarization vector field to treat the third 
harmonic generation. Roughly speaking, the study of the nonlinearity of the second 
order was aimed for presenting our method. The 'interesting' propagation equation 
system, (l5[ and ( fT6l l. show little interaction. In the third order, more interactions are 
considered, so that the physical effects are richer: the third harmonic generation, the 
optical Kerr-effect, cascade effects. To study Raman scattering, one has to generalize 
the result of this subsection to non-harmonic processes, as in done in [9|. 

A treatment similar to the one used to obtain the system (TBI) has to be done with 
now P = pW + P' 2 ) + P' 3 ). The propagation equation system we obtain is 
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M l p m (E p ) + (^) 2 (^LE g ,E p _ 9 j + £ LE 9 ,E r ,E p _ 9 _ r j) - -^J p . 

gez (<?,r)ez 2 £ ° 

(19) 

In agreement with the last paragraph of the last subsection, the interesting cases in 
the third order of nonlinearity are when the degree is also equal to three. We will thus 
consider only these cases. It is still assumed that the incident field is monochromatic, 
so that J p = JpSpiu. 

With the d\ definition In the d\ = 3 case, we have 

M i "(Ei)+o;f(2LE_ 1 ,E 2 J+3LE_ 1 ,E 1 ,E 1 j) =^-3 u (20a) 

M m (E 2 ) + (2wj) 2 |Ei,EiJ = 0, (20b) 
X! i m (E3) + (3 Wj ) 2 (2LE 1 ,E 2 J + [Ei.Ei.Eij) =0. (20c) 

This system presents several effects: the term |_E_i, E 2 J describes the counterac- 
tion of E 2 on Ei (as it appears with the second order of nonlinearity and the degree 

= 2), the pump beam is thus depleted; |_E_i, Ei, EiJ describes the optical Kerr- 
effect; lastly, the third harmonic is generated by two different processes: the cascade 
effect [Ei, E 2 J, and the third-order process [Ei,Ei,EiJ. 

With the doc definition In the doo — 3 case, we have 
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M in (Ei) (21a) 
+ w?(2LE_ 2 ,E 3 J +2LE_i,E 2 J 

+ 6|_E_ 3 , Ei, E 3 J + 3[E_i, E_i, E 3 J + 3[E_ 3 , E 2 , E 2 J 
+ 6LE_ 2 ,Ei,E 2 J +3LE_i,Ei,Eij) = ^-3 U 

M% n (E 2 ) (21b) 
+ (2 W/ ) 2 (2LE_i,E 3 J + LEi,E x J 
+ 6[E_ 3 , E 2 , E 3 J + 6[E_ 2 , Ei, E 3 J + 3[E_ 2 , E 2 , E 2 J 
+ 6LE_i,Ei,E 1 j) =0, 

X 3 m (E 3 ) (21c) 
+ (3w 7 ) 2 (2LEi,E 2 J 

+ 3LE_ 3 ,E 3 ,E 3 J +6LE_ 2 ,E 2 ,E 3 J +6LE_i,E!,E 3 J 
+ 3|E_ 1 ,E 2 ,E 2 J + L E ii E i) E ijl =0- 



The optical effects that this systems describes are the same than the ones present in 
the di — 3 degree: optical Kerr-effect, depletion of the pump beam, second and third 
harmonic generations. 

Once again, only energy considerations, presented in the following section, can 
justify that we keep terms like [E 2 ,E 2 ,E_ 3 J, that are, of course, several orders of 
magnitude smaller than |_Ei, Ei, E_iJ (see the equation (121 at "), 

This last example closes this section. The introduction of the degree allows to 
present systems whose solutions can be numerically investigated. Moreover, we see 
that, from the physical point of view (i.e., with the degree d\), the optical Kerr effect 
appears simultaneously to the depletion of the pump beam or the third harmonic gener- 
ation. This kind of statements is exactly the one we looked for; this is why we consider 
that, when restricting to harmonic processes with only one angular frequency in the 
incident field, we have fulfil our general aim. In (9) is exposed a method to generalize 
this when the scattered field oscillate at a nonharmonic frequencies or when the source 
oscillate at two frequencies. But this presents several drawbacks: first, this does not 
allow to treat phenomena described by a continuous set of frequencies, which are of 
course the most interesting ones from the physical point of view, and secondly, the 
systems obtained quickly become huge; solving them requires new hypotheses, so that 
it is like going back to the starting point: we have to impose constraints that depend on 
the particular situation we want to study. 

Remark 2 The optical Kerr effect 
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We said that the n-th harmonic generation is better studied in the n-th order of 
nonlinearity and the n-th degree; nevertheless, the third order of nonlinearity in the 
degree = 1 is also interesting. It consists of the single equation 

^"(EO+SwfLE-i.Ei.EiJ = =^Ji. (22) 

This is the optical Kerr effect. We repeat that we can argue that this system discards 
the term |_Ei,Ei,EiJ while keeping |_E_i, Ei, EiJ, thus 'neglecting' a priori that 
||E_i|| = ||Ei||. Nevertheless, in some cases, we can be concerned only to what 
happens in the fundamental frequency, and so we disregard the harmonic generation. 
This naturally leads to the study of this optical Kerr effect. We note that, dealing with 
this equation, we neglect the counter-reaction of the harmonics on the fundamental 
component, i.e., we neglect the term: 

2|_E_ 2 ,E 3 J +2[E_i,E 2 J + 6LE_ 3 ,Ei,E 3 J 

+3|_E_i, E_i, E3J + 3[E_3, E2, E2J + 6[E_2,Ei, E2J 

that contributes to P(wj). For non-centrosymmetric media, in which the second order 
polarization vector fields vanish, this approximation can be dangerous, since it is not 
always clear that 

2||[E_i,E 2 J|| «3||[E_i,Ei,EiJ||. 

However, for centrosymmetric media, the term | |_E_i, Ei, EiJ || must be the leading 
one among the terms that contributes to P^ 3 ' (loi); this justifies that the d x — 1 system 
is an approximation of the doo = 3 system. Hence, the way we obtained this equation 
shows that the optical Kerr effect appears naturally, and that it does make sense to 
study a harmonic nonlinear effect. 

We will go back, in the next section, to that simple system as a basic example of 
energy conservation. 

4 Electric energy considerations in nonlinear media 

4.1 The stochastic mean of the electric power density 

Before dealing with energy conservation in nonlinear optics, we have to give a rigorous 
definition of what we mean by energy transfer. Indeed, the components oscillating at 
the frequency of the incident field generate the other components, and so will loose 
some energy. The question is thus the following: when this harmonic generation occur 
with no loss of electric energy, that is, without Joule effect. Stated in an other way: 
when does the energy in the electric sources equal the energy in the scattered field? 

To this aim, we go back to the definition of the energy-momentum quadrivector, or, 
more precisely, its density, denoted by {cl, W). By definition, it is divergence-free, so 
that in a Cartesian coordinate frame, where 

V 4 - := (-cd x - cd y - cd z d t ); 
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we have 



-cV-(cX) + d t W = 0. 

From now on, we consider only the electromagnetic part (cX em , W em ) of (cl, W). 
We thus say that no transfer between electromagnetic energy-momentum and other 
form of energy-momentum occurs if 

-cV • (cl em ) + d t Wem = 0. 

For geometric reasons, and for the respect of the units of each field, the Poynting vector 
field, V := E x H, is defined to be the density of electromagnetic momentum cl em 
multiplied by the factor — c. The electromagnetic energy-momentum conservation is 
thus 

V • V + d t W em = 0. 

This relation allows to define the electromagnetic energy density W em . Indeed, Maxwell's 
equations give: 

d t W em = -V-V 

= (9 t B) • H + («9 t D) • E + E • J. 

We define the electric energy density W e by 

d t W e = (9 t D) • E. 

We write the electric vector field evaluated at the point s and the time t as (we 
need an absolute convergence of this series: indeed, in this case, the series converge 
uniformly with respect to the time variable and thus we will be allowed to change the 
integral operation with the sum operation) 

E(s,i) = ^E p (s)e-^ t . 

Once again, we will give all the details in the second order of nonlinearity, and then 
present more general formulae. So let us write 



D(t) = E E(t) + P(t) 

= e ]T (E„ + LE P J + £|E,,E p _,j)e- i **''«. 

The electric energy density is therefor^H 

1 Up to now, the dot product was on E 3 . We extend it on C 3 since we decompose the vectors in their 
Fourier components. We choose to take it linear in both variables, and so it will be symmetric - this will not 
bring problems since the nondegeneracy will never be exploited. 



20 



W e (t) = W e (0) 



(p,g)GZ*xZ P rSZ 

+ £ o Y -ipuiE-p ■ (E p + [E p J +^LE g ,Ep_,jJ)i. 

pez 9 ez 

The first term in this expression is a constant without any physical significance, the 
second one is oscillating and the last one diverges linearly with t. For large t (large with 
respect to the period associated to the fundamental angular frequency wj), only the last 
term is relevant. This leads to define the following quantity, which is the stochastic 
mean of the electric power density: 



<d t W e >:= lim i / dtd t W e 

T-s-oo T J Q 

= £ o Yl ( _ WUi^-v ■ ( e p + LEpJ + L E 9' e p-J )) • 



The term e (E p + [E p J + [Eg, E p _ 9 J), is the expansion of D p as £ E p + Pp + 

(2) 

Pp . We recall that, by construction, the answer of the medium is described by the po- 
larization vector field P; this is coherent with the fact that the sum — ipw/E-p • E p 

pGZ 

(where no susceptibility tensor appears) vanishes - we map p n> — p and use the Her- 
mitian symmetry of the duality product on C 3 . We thus have shown that, in vacuum, 
the stochastic mean of the electric power density e ^ — zpw/E-p • E p vanishes. 

pGZ 

Going back to an arbitrary order of nonlinearity, we have 



< 8 t W e > =< (dtP) • E > (23) 
= £ o ( _ ^lE-P ■ ( LEpJ + L E f E p-«J +•••)) 

pGZ qGZ 

neNpez 
= Y -ipuiE-p ■ Pp. 

pGZ 

We remark that this expression is a straightforward generalization of the stochastic 
mean of the electric power density in the linear regime. It is also clear that this expres- 
sion is real - taking the complex conjugate is equivalent to map p to — p (and q to — q 
in the second line), and thus does not change the sum. 
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4.2 A symmetry of the susceptibility tensors and lossless media 

We did not found a relation between the susceptibility tensors that exactly translate 
the property of a medium to be lossless. Nevertheless, we found one implication: if 
a certain symmetry between the susceptibility tensors is satisfied then the medium is 
lossless. This subsection is devoted to the presentation of this criterion. 

The first hypothesis that we will do is to consider that the electromagnetic power 
density vanishes in each order (we won't repeat 'the stochastic mean of each time, but 
it is always understood). We thus define 

pez 

The third equality in (|23l l, shows that 

< d t W e >=J2< d t W e n) > ■ 

Using the Hermitian symmetry of the Fourier component of the electric vector field 
and of the polarization vector field, it is straightforward to obtain the expression 

< d t W^ >= 2 W/ 3m{J] pE_ p • PW}, (24) 

4.2.1 In the first order 

In the first order, we have 

< d t W^ >: = e J2 -*PW/E_ P • LE p J 

= -ie^^P^-P ■ (x (1) {pui)-x (1)T {-puij)K p . 

pen 

Here we make the second assumption: we assume that each term of this sum vanishes. 
This leads to the conclusion that no electromagnetical energy is lost in the first order if 
the susceptibility tensor field (on M. 3 x R) x^ is Hermitian, i.e.: 

X W(p W/ ) = X (1)T (-p^). (25) 

The conclusion is not exactly the well-known fact from undergraduate lessons: if 
the medium is linear - so that the fields oscillate at only one angular frequency ojj -, it 
is lossless from the electrical point of view if and only if the susceptibility tensor field 

(on M 3 ) 

,Ldi) is Hermitian. Here, the medium is nonlinear - so that the fields 
oscillate at every angular frequency pojj for p in Z* -, and we showed that is lossless 
from the electrical point of view in the first order if the susceptibility tensor field 
is Hermitian. 
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4.2.2 In the second order 

The general criterion We now go to the second order: 

< d t W^ > = J2 -*PWjE- p • P {2) 

= £o X] -ipuiE-p ■ LE g ,Ep_gJ. 

(p,«)ez 2 

Some lines of computation show that 

O^r^q 

with 

:=r( X (2) (ro; / ,ga; / )-x (2)Tl2 (-(9 + r-V/.?^)) 
+ q( x {2) (rwj, gwj) - x (2) Tl3 (rwj, - (« + r)w 7 )) . 

where the partial transposition notation is defined by v a • Tl2 v&v c := Vb • v a v c 
for all triples (v a , Vb, v c ) of vectors in M 3 ; with indices, this gives (x' 2 ^ Tl2 )\ 1 i 2 = 
X^i/i 2 - In the same way we write v a • X^ Tl3v & v c := v c • x^ v i> v a and so 

(X^ 2 ^ Tl3 ) 1 j 1 j 2 = X^ 2 \ 2 ii- Also, d 2 (-, •) is the degeneracy function of a set of two 
elements: 

II, a = 6. 

For < 9tWe > to vanish, it is thus sufficient that each term of the sum vanishes 
(this is similar to what we did in the first order). We now introduce the third assump- 
tion: the factor of r and the one of q vanish independently. Thus, no transfer of electric 
energy to an other form of energy is possible if the tensor field x^ satisfies, for any 
(q, r) in 1? such that q > and < r < q, the relations: 

X (2) (««>/> = X (2)Tl2 (-{q + r)oJi,qui) 
X (2) (rwr,«wr) = x (2)Tl3 (r^i,-(q + r)oJi). 

A more symmetric way of writing these conditions appear if we define p := — (q + 
r) and we write explicitly, as an argument of the susceptibility tensor, the angular 
frequency for which x^ contribute in p( 2 ). Thus, from now on, any of the follow- 
ing expression has the same meaning: x^ n \ s 'i w o! • • • x'"'(s,wi,'" 
x'™-* (wo; Wi, • ■ • , u> n ), x^™H w i; ' ' ' i w n)> an d when x*-"-* is evaluated onn + 1 angular 
frequencies ojo, • • • , w n > it is understood that cjq = c*>i + • • • + cj„. With this notation, 
a medium is lossless from the electrical point of view in the second order if 
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X (2) (-pwr;r-wr,<7wr) = x (2) Tl2 {-rui;pui, qui) = X {2) Tl3 (-<?<*>/; rui,pui). 

(26) 

So no electric energy loss is guaranteed if the symmetry that consists in permuting the 
indices together with the frequency variables holds for the x ■ This result appears 
in the book [10|, but these authors do not consider the cascade effects that appear in 
harmonic generations. Moreover, the assumptions needed to derive these results are 
not explicitly given. 

We note that (contrary to what is affirmed in |2|, p. 34) we can have a lossless 
medium with complex- valued susceptibility tensor fields. This has been checked nu- 
merically by the authors, with simulations similar to the ones given in [3 |. 



Kleinman's relations If the medium is instantaneous, then the response function 
iZW is of the kind A 5(U) where A is a function that depends only on the 

ie{l,-" ,n} 

space variables. Then the susceptibility tensor fields do not depend on the angular 
frequencies: 

xtl ■= x {n) (Piui, ■ ■ ■ , Pn Ui) V(pi, • • • , Pn ) e Z n . 
Thus, in the first order, the polarization is 



/CO 
-oo 

/oo 
dwiE(wi)e 
-oo 



-iwit 



insi 

= £o xi 1 l ) st E(i). 

Since Xinsv wnen contracted with a real vector, gives a real vector, it is real. 
In the same way, in the second order, we have 



/oo poo 
dwi / ^ 2X ^ st E( Wl )E( W2 )e-^ +I 
-oo J — OO 

/OO pOO 
rf Wl E( Wl )e- lWlt / du 2 E(u 2 )e 
-oo J — oo 



and Xinst 1S also real - 

Reporting the criteria ( T25l ) and d26l i in this context, we deduce that an instantaneous 
nonlinear medium is lossless if the susceptibility tensor fields that characterize it satisfy 

f in the first order : X flt = xfll 

\ in the second order : x^L = XinJ™ = xfl^ 13 ■ 
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Because the transposition (12) generate the symmetric group S2, the transpositions 
(12) and (13) generate £3, etc., where <S„+i acts as permutation of the components of 
X^ n \ we conclude that completely symmetric and frequency independent susceptibility 
tensor characterize lossless media. This is the famous result of D.A.Kleinman ( |[TTl0 . 

In the lower degrees In the second order 1 < n < 2 and the degree doo = 2, we 
have 

P^= £ oLE 2 J, 
Pf } =2eoLE_ 1 ,E 2 J, 
=e LEi,E 1 J. 

Hence, an explicit expansion of (l24l i gives 

< d t W e >= 2eo W/ 3m{E_ t • ( LE, J + 2 LE_, , E 2 J ) + 2E_ 2 • ( LE 2 J + [E^EjJ)}. 

We first treat the first order of nonlinearity. So let us assume that x' 1 ) (s, wj) is Hermi- 
tian for all s in M 3 , in the sense that its transpose conjugate is equal to itself; we then 
have 

E_! • [ExJ : = E_ 1 - X (1) (^/)Ei 
= E 1 - X (1)T (^/)E_ 1 

= e _ 7-x (1) (^/)eT 

= E_ 1 - X (1) (^/)Ei 
=:E_i- |EiJ, 

We deduce from this that 3m{E_j • [EjJ} vanishes: this is the condition (f25t on 
electrical energy conservation in the first order. A similar computation shows that if 
X^( s > is Hermitian for all s in M. 3 , then 3m{E_j • [EgJ } vanishes. 

Considering now the second order of nonlinearity, we suppose that the second order 
susceptibility tensor field satisfies 

X (2) (~pu>i;ruj I ,quj I ) = x i2)Tl3 {-qur,ruj I ,pLL> I ) (27) 

for all [q, r) in Z 2 - we will use only x^ Tl3 (^i; —u>i, 2u>i) — x^ (— 2w/; — — oji). 
Then we have 
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< d t Wf> > = 4e w / Sm{E_ J • [E^EgJ + E_ 2 • LEi,EjJ} 

= 4e w / 3m{E_ J • X (2) (-^i, 2w/)E_ j E 2 + E_ 2 • X {2) (^i, ®i } 

= 4e w/3^{E 2 •x (2)Tj3 (-^/^^/)E-jE_ j + E_ s • x (2) (w/, w/)E J E J } 
= 4e w/3m{E 2 -^'(-w/^wjJE^E^ + E_ 2 • x (2) (w/, w/JEj Ej } 

= 8e ui3m{«e{E s ■ X [2 \-^i, -w/)E_ J E_ i }} 
= 0, 

where the last but one line is obtained by the Hermitian symmetry of the electric vector 
field and of the susceptibility tensor fields. We thus see that the criterion d27| i is a 
sufficient condition for the stochastic mean of the electrical energy density at the second 
order to vanish. 

This is why we argued, in the subsection !3.2l that though the degree d\ is the most 
satisfying one from the physical point of view, even if it leads to consider Ei as a tank, 
the degree doo allows to consider lossless media. This advantage can be a good test to 
check the coherence of numerical simulations. 

Lastly, we note that d27l i is weaker than d26b ; this is due to the degeneracy in the 
angular frequencies in the order two and degree two. The reader can check that the full 
criterion d26i i is necessary if we consider the order two in the degree three. 

4.2.3 In the third order 

The general criterion In the third order, we have 



then a medium is lossless from the electrical point of view in the third order. 

In the lower degrees As exposed in the remark|2] the system of propagation equa- 
tions ( fT9b in the degree d^ = 1 reduces to the single equation that describes the optical 
Kerr effect: 



< d t W^ >:= e -*8f/E_, ■ LEr,E s ,E g _ r _ s J. 



We can show (|9|) that if the following conditions are satisfied 



< x^i-p^iiquJi, ru)i,su>i) = x (3) Tl3 {-ruir, qwi,pu>i, sui) 
^(3) (—pun qu^rwi, aujj) = x (3) Tli {-suji]quji,ruji.puji) 



Mi in (Ei) + 3a;|LE-i,Ei,EiJ 



— l UJl 



Ji. 



((EH) 



Since 




£o LE-i, Ei, EiJ 
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the formula d24l > gives 



pSN 

= 6e w/3m{E_j • LE_j,Ei,EjJ} 

Assume furthermore that the medium is invariant along one axis (say (O, z)), and 
that if a beam linearly polarized along (O, z) impinges on this medium, then the scat- 
tered light is also polarized along (O, z). Then, the only components of the suscep- 
tibility tensors that matter are x^ Z * an d X^ Z zzz- Let us denote Ei(a;, y, z) by 
u\{x,y)z. 

In this case, we have 

< dtW^ > K eir = 6e w/3m{M_ J x (3)z zzz (-w/,w/,w/)w-jwjttj} 
= &e () uj I '^m{x {3)z zzz (-w/, w/, \ 4 . 

In other words, with all these assumption, there are no electrical energy loss in the 
third order if and only if the susceptibility tensor x^ is real-valued. 

4.3 In the n-th order 

The general proof that a sufficient condition for a medium to be lossless in the n-th 
order can be found in J5J- The result is that 

pel. 

= -ie ux ^2 + ^P'i) e -(pi+---+p„) ' L e pi; ' ■ ' ; E pJ- 

vanishes if 

X^ Tli+1 (-PjUr,PaUi,--- ,Pj-iui,Pj+iui, ■ ■ ■ ,p n ui) = 

x in) {-p un;piui, ■ ■ ■ ,p n uji) Vj e {1, • • • ,n} 

where po := — (pi + ■ ■ ■ + p n ). We can check that the criteria given for the lowest 
order fulfils this set of conditions. 

Once again, we repeat that we did not find a necessary and sufficient criterion: we 
supposed that each order vanishes independently, that the term < <9 t ( We™" 1 > vanishes 
for each set (-po,pi, ■ ■ ■ ,p n ), withp = -(pH hp n ). where E po , E Pl , • ••, E Pn 

(n) 

appear in E_ p • P. p , and finally that, within each such set, the factors of p\, etc., p n 
vanish separately. These factors are of course the terms 

X ( " )Tli+I (~Pj w i;PoWr,-- - ,pj-iu I ,pj + - i ,u)i,--- ,p n ui) 
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-x (n) (-Pow/;piw/, • • • ,PnUl), 

justifying the general criterion. 

This general result appears in iflOl ; we hope that the physical meaning of these 
equations, as well as the assumptions we were led to carry out, have been clarified. 

5 Conclusion and outlook 

This paper explores a new route to obtain the propagation equation systems governing 
nonlinear interactions of light in matter. The aim is to obtain equations that do not 
depend on the particular effect, the particular involved material or the particular situa- 
tion we want to study. This aim is partly fulfilled, especially when treating harmonic 
generation. For this, the introduction of the degree is of major importance. To our 
knowledge, no similar notion has been introduced in the literature. 

The paper [3| is devoted to numerical results; the system d2TT > is chosen for its 
interesting effects (second and third harmonic generation, optical Kerr effect, depletion 
of the pump beam) and because it allows to make a power balance. Of course, only 
an experimental test could invalidate the theoretical as well as the numerical models. 
Though these experiments should not be difficult to carry out, we have to confess, 
without being to much ironic, that the result is highly uncertain. Indeed, we fear to fall 
into one of the two possibilities: either the incident power is too low for the nonlinear 
effects to be measurable, or, particularly because our incident beam is monochromatic, 
the medium will not be insensitive to light, so that thermal or mechanical effects should 
be included. 

To conclude, the method exposed in this paper has to be generalized in several 
directions: to match experimental data, we have to deal with non-stationary media, to 
include non-harmonic generation (in order to study subharmonic generation, Raman 
scattering, etc.), and also to generalize our notion of the degree to the cases where the 
sources of the electromagnetic field oscillate at several frequencies (in order to study 
sum-harmonic generation, four-wave mixing, etc.), or even with a continuous spectrum 
(in order to deal with laser pulses). 
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